home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dlassq.z / dlassq
Encoding:
Text File  |  2002-10-03  |  3.4 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAASSSSSSSSQQQQ((((3333SSSS))))                                                          DDDDLLLLAAAASSSSSSSSQQQQ((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLASSQ - return the values scl and smsq such that  ( scl**2 )*smsq = x( 1
  10.      )**2 +...+ x( n )**2 + ( scale**2 )*sumsq,
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DLASSQ( N, X, INCX, SCALE, SUMSQ )
  14.  
  15.          INTEGER        INCX, N
  16.  
  17.          DOUBLE         PRECISION SCALE, SUMSQ
  18.  
  19.          DOUBLE         PRECISION X( * )
  20.  
  21. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  22.      These routines are part of the SCSL Scientific Library and can be loaded
  23.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  24.      directs the linker to use the multi-processor version of the library.
  25.  
  26.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  27.      4 bytes (32 bits). Another version of SCSL is available in which integers
  28.      are 8 bytes (64 bits).  This version allows the user access to larger
  29.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  30.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  31.      only one of the two versions; 4-byte integer and 8-byte integer library
  32.      calls cannot be mixed.
  33.  
  34. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  35.      DLASSQ returns the values scl and smsq such that ( scl**2 )*smsq = x( 1
  36.      )**2 +...+ x( n )**2 + ( scale**2 )*sumsq, where  x( i ) = X( 1 + ( i - 1
  37.      )*INCX ). The value of  sumsq  is assumed to be non-negative and  scl
  38.      returns the value
  39.  
  40.         scl = max( scale, abs( x( i ) ) ).
  41.  
  42.      scale and sumsq must be supplied in SCALE and SUMSQ and
  43.      scl and smsq are overwritten on SCALE and SUMSQ respectively.
  44.  
  45.      The routine makes only one pass through the vector x.
  46.  
  47.  
  48. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  49.      N       (input) INTEGER
  50.              The number of elements to be used from the vector X.
  51.  
  52.      X       (input) DOUBLE PRECISION array, dimension (N)
  53.              The vector for which a scaled sum of squares is computed.  x( i )
  54.              = X( 1 + ( i - 1 )*INCX ), 1 <= i <= n.
  55.  
  56.      INCX    (input) INTEGER
  57.              The increment between successive values of the vector X.  INCX >
  58.              0.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAASSSSSSSSQQQQ((((3333SSSS))))                                                          DDDDLLLLAAAASSSSSSSSQQQQ((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      SCALE   (input/output) DOUBLE PRECISION
  75.              On entry, the value  scale  in the equation above.  On exit,
  76.              SCALE is overwritten with  scl , the scaling factor for the sum
  77.              of squares.
  78.  
  79.      SUMSQ   (input/output) DOUBLE PRECISION
  80.              On entry, the value  sumsq  in the equation above.  On exit,
  81.              SUMSQ is overwritten with  smsq , the basic sum of squares from
  82.              which  scl  has been factored out.
  83.  
  84. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  85.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  86.  
  87.      This man page is available only online.
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.